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Quantum Phase Transitions in Two-Dimensional Spin Systems with 
Ladder, Plaquette and Mixed-Spin Structures 
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Quantum phase transitions between the magnetically ordered and disordered states are stud- 
ied for the two-dimensional antiferromagnetic quantum spin systems with ladder, plaquette, 
and mixed-spin structures. Starting with properly chosen singlet-cluster configurations, we per- 
form the series expansion for the staggered magnetic susceptibility. The phase boundary is 
determined by applying the Dlog and biased Pade approximants to the staggered susceptibility 
thus obtained. The resulting phase diagram allows us to discuss the quantum phase transitions 
quantitatively, which agrees fairly well with the quantum Monte Carlo results for several cases 
previously studied. 
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§1. Introduction 

Recently low-dimensional spin systems with the gap 
for the excitation spectrum have been the subjects of 
considerable interest. The spin gap is generated by var- 
ious mechanisms, which may be sensitive to the lattice 
structure, the competing interactions, the topological na- 
ture of spins, etc. A typical example is the spin plaquette 
system such as CaV409jJ'aLl'Q'ElQlZl'EP which may be de- 
scribed by the two-dimensional (2D) Heisenberg model 
on a Ij^ 5.^ depleted square lattice or the meta-plaquette 
modelBEP The disordered ground state with the spin gap 
observed in these systems may result from the singlet- 
spin configuration formed in each plaquette. Another 
prototypical example is the coupled spin-laddejc-SistemEP 
realized in the coipcounds such as SrCu2 03ll3'tlP and 
(Sr, Ca)i4Cu2404i.tHr For such ladder systems with even 
number of legs, the spin gap is generated to stabilize the 
disordered ground state. Interestingly enough, the in- 
troduction of non-magnetic impurities into this system 
induces iba^hase transition to the magnetically ordered 
state.EJ'M'H Furthermore, to the family of the spin gap 
systems we may add the mixed-spin system, in which the 
topological nature of spins as well as the lattice structure 
play an essential role to produce the spin gap. For in- 
stance, in one-dimensional (ID) systems, a mixed-spin 
chain with the alternating array of two kinds of spins 
has stiniiilatod intensive experimentalo^ and theoreti- 
cal0-0B0@ studies. It is known that the ID mixed- 
spin chain realizes either the ferrimagnetic or the singlet 
ground state depending nn bnw wp arrange diiferent type 
of spins on the ID latticeH'B0 

Various spin gap systems mentioned above provide 
us with an interesting research area of quantum spin 
systems. Among others, the quantum phase transition 
from the spin gap phase to the magnetically ordered 
phase is one of the most interesting issues. Motivated 
by these hot topics we investigated in the previous pa- 



peicj) the competition between the magnetically ordered 
and disordered states in the 2D quantum spin system 
which includes both of the bond- and spin-alternations. 
By means of the non-linear signia model and the modi- 
fied spin wave approach, we evaluated the spin gap and 
the spontaneous staggered magnetization to discuss the 
quantum phase transitions. Although this study enabled 
us to qualitatively describe how the spin gap phase is 
driven to the antiferromagnetic one, the resulting phase 
diagram was not sufficient enough to give quantitative 
discussions. 

The purpose of the present work is to study more 
quantitatively the quantum phase transitions for the 
2D spin systems with ladder, plaquette, and mixed-spin 
structures. To this_epd we employ the series expansion 
method developedrHPT and extensiiieht used by many 
groups far ladders,E3) 2D sjjstems,^^^^ Kondo lattice 
modelsJSa' bilayer systems,E3) etc. In this approach, 
starting with properly chosen singlet-spin clusters, we in- 
troduce the couplings among the clusters perturbatively 
to carry out the series expansion of the physical quanti- 
ties such as the staggered susceptibility. Then the critical 
point for the phase transition is determined by applying 
the Pade approximantsEa^ to the physical quantities thus 
obtained. 

The paper is organized as follows. In §2, we first in- 
troduce the Hamiltonian for 2D spin systems and outline 
how to apply the series expansion techniques to our sys- 
tems. In §3, by performing the series expansion for the 
staggered susceptibility and then employing the Pade ap- 
proximants, we obtain the phase diagram to discuss the 
competition between the disordered and ordered states 
for three kinds of 2D quantum spin systems mentioned 
above. Brief summary is given in the last section. 

§2. Cluster Expansions 

We consider a 2D antiferromagnetic quantum spin 
system, which is described by the following generalized 
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Heisenberg Hamiltonian 
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where Sj is the spin operator at the j-th site on the 
square lattice and Ji, J2, and J3 are the antiferromag- 
netic exchange coupUngs (Ji, J2, J3 > 0). Note that the 
spin Sj is allowed to take different values at each cite, 
which enables us to treat the systems with mixecLspins. 

In order to apply cluster expansion techniques tJ^ the 
Hamiltonian is divided into two parts. Since we start 
with a strong-coupling singlet state, we take the unper- 
turbed Hamiltonian Hq as an assembly of independent 
singlet-spin clusters formed by the coupling Ji (the cor- 
responding bonds {i,j) are sepecified by Di). The inter- 
actions among independent clusters are then taken into 
account by series expansions in the perturbed term Hi. 
We introduce two types of perturbations, for which the 
corresponding sets of bonds {i,j) are denoted by D2 and 
D3. The detail of Di,D2, and D3 will be given for each 
case in the following sections. We henceforth assume 
that A(= J2/J1) < 1 and aA(EE J3/J1) < 1(0 < a < 1). 
We will see that this parameter regime indeed includes 
physically interesting cases. 

In the following, we discuss the quantum phase tran- 
sitions for the 2D antiferromagnetic spin systems with 
ladder, plaquette, and mixed-spin structures. To this 
end, we introduce three kinds of singlet-cluster configu- 
rations, i.e., the dimer singlet, the plaquette singlet, and 
the mixed-spin-cluster singlet. The last one, which is 
unique for our systems, is composed of a specific singlet- 
spin configuration, e.g. 1/2 o 1 o 1/2 (see Figs. 8 and 9). 
The introduction of the mixed-spin-cluster allows us to 
perform a systematic series expansion for 2D mixed spin 
systems. Starting with the above spin singlet states, we 
can carry out the series expansion with respect to A as 
well as aX. For the ladder, the plaquette, and the mixed- 
spin structures, we refer to the corresponding expansions 
as the dimer, the plaquette, and the mixed spin-cluster 
expansions, respectively. We naturally expect that the 
introduction of Hi perturbs the singlet ground state with 
the excitation gap and gradually enhances the antiferro- 
magnetic spin correlations, finally giving rise to the long- 
range magnetic order. This will be shown to be indeed 
the case for our systems. 

We calculate the staggered spin susceptibility to deter- 
mine the critical point for the transition. We thus add 
the following Zeeman term as a perturbation. 
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where h is the staggered magnetic field and A (B) de- 
notes one of the two sublattices. We estimate the ground 
state energy E{h) of the total Hamiltonian H + HgT up 
to the second order in h, and then obtain the magnetic 



susceptibihty x for the staggered field. 
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This quantity is expanded as a power series in A as 



x-E«"^"- 
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(2.6) 



We calculate the susceptibility up to the eighth order 
in A for the dimer expansion and the fourth order for 
both of the plaquette and the mixed-spin cluster expan- 
sions. We should recall here that the phase boundary 
between the magnetically ordered and disordered states 
is given by the critical line on which the staggered sus- 
ceptibility is divergent. Therefore a further approximate 
procedure is necessary to deduce the sensible singularity 
by the asymptotic analysis of the power-series expansiou. 
For this purpose, we make use of Fade approximantsESP 
for the susceptibility obtained up to the finite order in 
A. Besides ordinary Dlog Fade approximants, we also 
employ biased Fade approximants,ES' for which we as- 
sume that the phase transition in our 2D quantum spin 
models should belong to the^universality class of the 3D 
classical Heisenberg model. Ej^ Then the critical value of 
Ac is determined by the facpula X ^ (A — Ac) •" with the 
known exponent 7 = 1.40 We shall see below that the 
biased method provides a fairly good approximation for 
Ac in some cases, and in general is useful to check how 
well our Fade approximants work. 

§3. Quantum Phase Transitions 

In this section, we discuss the quantum phase tran- 
sitions by applying the series expansion to our models. 
We have two parameters in the perturbed term, A and 
a, so that we can observe in different ways how the 
2D antiferromagnetic correlations develop. To confirm 
the validity of our series expansions and Fade approxi- 
mants, we first study the 2D system with ladder struc- 
ture which was already studied extensively by various 
methods. Our results are compared with those of the 
quantum Maiite|i^arlo (QMC) simulation in rather good 
agreement.EilBEd) We then move to the plaquette spin 
systems and the mixed spin systems, and argue how the 
plaquette-singlet state and the mixed-spin singlet state 
are driven to the 2D magnetically ordered state. 



3.1 Ladder- structure systems 

Let us start with a 2D spin system with the ladder 
structure, which is shown schematically in Fig. 1. In 
this figure, the circle represents the s — 1/2 spin sit- 
ting on the square lattice. Note that this 2D model 
can be constructed from isolated dimers in two ways de- 
pending on how we may introduce the couplings among 
dimers perturbatively: we may refer to the system corre- 
spondingly as the coupled 2-leg ladders and the coupled 
dimer chains. In the former case, the bold, the thin, 
and the dashed lines in Fig. 1 represent the coupling 
constants 1(= Ji),Al(= J2), and aLAL(= J3), respec- 
tively. By tuning the value of Ol, we naturally inter- 
polate the independent 2-leg ladders (c^l = 0) and the 
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Table I. Series coefficients an for the dimer expansion of the stag 
gered susceptibihty per site x for the coupled-ladder system. 
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Fig. 1. 2D spin system with the ladder structure. The circle 
represents the spin s = 1/2. We refer to the system as the 
coupled-ladder (coupled-dimer-chain) system, when the bold, 
the thin, and the dashed lines represent the coupling constants 
1, AL(aD'^D)i and cti^Xi^^Xi)), respectively. 



2D systems. On the other hand, in the latter case, by 
taking 1(= Ji),Ad(= J3), and aDAD(= J2), we can in- 
vestigate how the independent dimer chains {a^, — 0) 
are combined to make the 2D systems. We distinguish 
these two constructions which may be complementary 
to each other because the available parameter regime is 
restricted in our series expansion approach. 
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Fig. 2. Phase diagram for the coupled-ladder system in Fig. 1. 
The solid (dashed) line indicates the phase boundary obtained 
by the biased [3/4] (Dlog [3/4]) Pade approximantsj-^iber^lcd 
circles represent the results of the QMC simulationsJ2irE3'E3' 



Let US first regard the 2D system as the coupled 2-leg 
ladders as shown in Fig. 1. The starting Hamiltonian if 
has the ground state with the spin gap due to the dimer 
singlets. We calculate the staggered susceptibility x by 
means of the dimer expansion up to eighth order in Al for 
various values of aL . The resulting power series for some 
particular values of aL are presented in Table I. Using 
Pade approximants, we get the phase diagram shown in 
Fig. 2. In this figure, the solid (dashed) line represents 
the phase boundary obtained by the biased [3/4] (Dlog 
[3/4]) Fade approximants. When ol = 0, the system 
is reduced to the isolated 2-leg ladders with the interleg 
(intraleg) coupling constant 1 (Al), which Js known to 
have disordered ground state with spin gap.B* Increasing 
the parameter q;l with a fixed Al , the antiferromagnetic 
correlation grows up, and eventually the phase transition 
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11.400106 


37.401693 
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1.3427618 
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0.32560333 
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30.989600 
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3.6797174 


43.252600 
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to the antiferromagnetically ordered state occurs. For 
instance, if we determine the phase boundary by means 
of biased [3/4] Fade appoximants, the critical value is 
given by aL — 0.26 for Al — 1. 

We find that near q;l ~ 0.8 the phase boundary deter- 
mined by the Dlog Fade approximants exhibits a patho- 
logical behavior, namely it branches out into upper and 
lower lines. It may be obvious that these lines may not be 
physically sensible. It is known that this type of pathol- 
ogy occasionally appears in Fade approximants. If we 
discard these spurious parts in critical lines, we then find 
that the results of two Fade approximants show the com- 
mon behavior and are both in good agreement with those 
of the QIVIC simulations.E^'L3'E3 Also, our results repro- 
duce the critical value Al = 0.54 fpx q;l = 1 which was 
previously obtained by Singh et a/.cil' We should notice, 
however, that the biased Fade approximants may not al- 
ways give more accurate results than the ordinary Dlog 
Fade approximants. We must carefully determine the 
phase boundary after trying various Fade approximants, 
as will be momentarily shown below. 
We next regard the present system as the coupled dimer 
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Fig. 3. Phase diagram for the coupled-dimer-chain system in Fig. 
1. The solid (dotted) line indicates the phase boundary obtained 
by the Dlog [4/3] (biased [4/3]) Pade approxjiBri|a«ts. The filled 
circles represent the QMC simulation results.Eil'Ej' 



chains, for which the corresponding couplings are defined 
in Fig. 1. Repeating a similar calculation as in the previ- 
ous case, we arrive at the phase diagram shown in Fig. 3. 
We list the obtained series for some values of a in Table 
II. In this figure, the solid and the dotted lines represent 
the phase boundaries obtained by the Dlog [4/3] and 
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Table II. Series coefficients a„ for the dimer expansion of the 
staggered susceptibihty per site x fo'' the coupled-dimer-chain 
system. 

n qd = 0.0 od = 0.2 ao = 0.5 qd = 1-0 
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Fig. 4. 2D spin system with the plaquette structure. The circle 
represents the spin s = 1/2. The bold, the thin, and the dashed 
lines represent the coupling constants 1, A, and cxX. 



the biased [4/3] Pade approximants. In the latter anal- 
ysis, we have two critical lines, one of which (labeled by 
crosses) exhibits a pathological behavior quite different 
from that of the Dlog Pade analysis. If we discard this 
line, we find that the results of two Pade approximants 
show the physically sensible bcha\doiL^nd are consistent 
with those of the QMC simulationS'Ej) (the filled circles 
in Fig. 3). 

Finally we wish to note that at the point (q;d,Ad) ~ 
(0, 1), our system just lies on the critical line which sep- 
arates the ordered and disordered states, as seen from 
Fig. 3. Since the system in this case is reduced to the 
independent isotropic spin chains, our numerical results 
reproduce the well-known fact that the ground state of 
the spin-1/2 Heisenberg chain is in a critical spin liquid 
phase, so-called Tomonaga-Luttinger liquid nliase, with 
neither the spin gap nor the long-range order .EiP We also 
note that the properties around this quantum ^critical 
point have been already studied by Affleck et al^B 

The above studies on the ladder-structure systems 
seem to give rather satisfactory results, which encour- 
age us to apply a similar series expansion approach to 
the analysis of other quantum spin systems. 

3.2 Plaquette- structure systems 

Let us now turn to the plaquette spin systems. We 
here consider two kinds of the models. First, we are con- 
cerned with the system shown in Fig. 4. The starting 
Hamiltonian /Jq describes a sum of the independent pla- 
quettes, whose ground state is spin singlet with the exci- 
tation gap. Introducing the perturbed part Hi may in- 
duce the competition between the plaquette-singlct cor- 
relation and the antiferromagnetic one. By this construc- 
tion of the system, we naturally interpolate the isolated 
plaquettes A = (or isolated ladders a = 0) and the nor- 
mal square lattice. The coefficients are calculated up to 
the fourth order in A by means of the plaquette expan- 
sion and tabulated in Table III for some particular values 
of a. Using the biased [1/2] Pade approximants, we ob- 
tain the phase diagram shown in Fig. 5. In this phase 
diagram the system on the a = axis is composed of the 
assembly of the independent 2-leg ladders, which belongs 
to the disordered spin-gap phase. Away from this axis, 
the 2D antiferromagnetic correlation grows up and the 
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Fig. 5. Phase diagram for the plaquette system in Fig. 4. The 
solid line indicates the phase boundary obtained by the biased 
[1/2] Pade approximants. 



Table III. Series coefficients a„ for the plaquette expansion of the 
staggered susceptibility per site x for the 2D spin system with 
the plaquette structure. 



: 0.0 



:0.2 



0.5 



1.0 



1.3333333 1.3333333 1.3333333 1.3333333 

1 1.7777778 2.1333333 2.6666667 3.5555556 

2 1.6009195 2.6131044 4.3715197 7.9425797 

3 1.1215579 2.9136374 6.8339620 17.102341 

4 0.5614964 3.0347609 10.187867 35.146612 



quantum phase transition to the ordered state occurs at 
the critical value ac- When A = 1, we find ac = 0.22. 
We note that this plaquette system with A = 1 is equiv- 
alent to the coupled 2-leg ladders with Al = 1 discussed 
in the previous subsection, where we have obtained the 
slightly different critical value ol = 0.26. It is also to be 
noticed that the mean-field theory by the bond-operator 
representation and the QMC simulation also-yield the 
corresponding values ac = 0.25E3' and 0.32E3) for the 
system with A = 1. These results imply that it may be 
necessary to carry out higher order cluster expansions 
both for dimer and plaquette systems to obtain more ac- 
curate critical values of a in the region A J^ 1 (This is 
indeed seen from Fig. 2 for the ladder system). Finally 
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we wish to mention that for the special case (a = 1), 
where each plaquette is connected with nearest neighbor 
ones via the single coupling constant A, similar results 
have been reported by Fukumoto and Oguchi. 




Fig. 6. 2D spin system composed of the plaquette chains. The 
circle represents the spin s = 1/2. The bold, the thin, and the 
dashed lines indicate the coupling constants 

1, A, and a A. 
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Fig. 7. Phase diagram for the plaquette system in Fig. 6. The 
solid line represents the phase boundary obtained by the biased 
[2/2] Pade approximants. 



Let us next consider the system shown in Fig. 6, 
which mayJ^e considered to be made out of plaquette 
chains,QE3'0 because this model, for a = with fi- 
nite A, is reduced to the isolated chains with plaquette 
structures. An interesting point is that this system is 
toppJjQgjcally equivalent to the 1/5 depleted square lat- 
ticeBQ'© which is relevant to the study for the compound 
CaV409. The phase diagram shown in Fig. 7 is obtained 
by the fourth order calculation and the biased [2/2] Pade 
approximants. The resulting series for some values of a 
are listed in Table IV. Let us first look at the results with 
the value of A being fixed. By increasing the value of 
a from zero, we observe the evolution of the plaquette 
chain to the depleted square lattice. For example, in 
the case of A = 1, we get the critical value ac = 0.82. 
On the other hand, when we fix a = 1, we can see how 
the isolated plaquettes are uniformly coupled to form 
the 1/5 depleted square lattice, for which we obtain the 
critical value Ac = 0.91. This value is in goodj^gree- 
ment with the result already oblained by QMCa) and 
also by the plaquette expansion.EP As clear from these 



Table IV. Series coefficients a„ for the plaquette expansion of 
the staggered susceptibility per site x for the 2D spin system 
composed of the plaquette chains. 
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2.1449010 


3 


0.20889382 


0.47565756 
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2.6268926 
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0.073161406 


0.29653125 


0.92711496 


3.1394401 



results, the plaquette structure is quite essential, but is 
not sufficient to produce the spin gap for the isotropic 
1/5 depleted square lattice (A = a = 1), for which we 
have the antiferromagnetic long-range order. So, it is 
necessary to introduce the dimer structure Qr_frustrating 
couplings to have the spin gap in this case.oQ-' 

3.3 Mixed-spin systems 

So far, we have restricted our discussions to the 2D 
spin-1/2 systems, for which the lattice structure plays 
an important role to generate the spin gap. As men- 
tioned in the introduction, the mixed-spin systems with 
periodic array of different spins also have attracted much 
attention both experimentally and theoretically. In these 
systems, not only the lattice structure but also the topo- 
logical nature of spins become imnartant for the gap for- 
mation. In the previous papers,t£lL3' we have investi- 
gated how the ID mixed spin chains are coupled to form 
the 2D mixed spin system by means of the non-linear a 
model and the modified spin wave analysis. Although 
the quantum phase transitions between magnetically or- 
dered and disordered phases were described at least qual- 
itatively by the above approaches, the results obtained 
turned out to be far from quantitative discussions. The 
purpose in this subsection is to quantitatively explore 
the phase transitions in 2D mixed spin systems by using 
the series expansion method. To be specific, we wish to 
deal with two typical systems composed of s = \/2 and 
s = 1, which are shown in Figs. 8 and 9. 

Let us first consider the system shown in Fig. 8, which 
we refer to as the " columnar mixed-spin system" , for 
which the mixed spin chains arc stacked uniformly in a 
vertical direction. In this figure, the small circle (large 
filled circle) represents s — l/2(s = 1) and the bold, the 
thin, and the dashed lines indicate the coupling constants 
1(= Ji),A(= J2), and a\{= J3), respectively. When 
A = 0{Hi — 0), the wave function of the ground state is 
the direct product of mixed-spin-cluster singlets given by 
the spin arrangement of 1/2 o 1 o 1/2. This is our starting 
configuration for the cluster expansion. The phase tran- 
sition to the ordered state may be anticipated when the 
values of a and A are increased. We perform the mixed- 
spin cluster expansion up to the fourth order and list the 
series coefficients for some particular values of a in Ta- 
ble V. Employing the biased [2/1] Fade approximants for 
the staggered susceptibility, we end up with the phase 
diagram shown in Fig. 10 (solid line). We find that the 
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Fig. 8. 2D spin system with the columnar mixed-spin structure. 
The small (large filled) circle represents s = 1/2 {s = 1). The 
bold, the thin, and the dashed lines indicate the coupling con- 
stants 

1, A, and aX. 




Fig. 9. 2D spin system with the diagonal mixed-spin structure. 
The small (large filled) circle represents s = 1/2 (s = 1). The 
bold, the thin, and the dashed lines indicate the coupling con- 
stants 

1, A, and aA. 
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Fig. 10. The phase diagram for the systems in Fig. 8 (the solid 
line) and Fig. 9 (the dashed line). The critical lines are deter- 
mined by the biased [2/1] Fade approximants. 



mixed spin chains {a — 0) are stacked diagonally. The 
definition of the coupling constants is the same as that 
in Fig. 8. The mixed-spin cluster expansion up to the 
fourth order with the biased [2/1] Fade approximants 
yields the phase diagram shown in Fig. 10 (the dashed 
line) . We tabulates the resulting series for some values of 
a in Table VLFor a = 0, the system is correctly reduced 

Table VI. Series coefficients a„ for the mixed-spin cluster expan- 
sion of the staggered susceptibility per site x foi* the diagonal 
mixed-spin system. 



0.0 



: 0.2 



:0.5 



1.0 



1.7777778 1.7777778 1.7777778 1.7777778 

1 1.1851852 2.3703704 4.1481481 7.1111111 

2 0.63981053 2.8229887 8.3587481 23.614326 
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Table V. Series coefficients a„ for the mixed-spin cluster expan- 
sion of the staggered susceptibility per site x for the columnar 
mixed-spin system. 
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4.0253708 
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88.590346 



0.097548541 4.7431208 38.497553 259.24615 



system at the point {a, A) = (0, 1), which consists of in- 
dependent mixed-spin chains, is in a disordered phase 
with spin gap. This result correctly reproduces the fact 
that the low-energy excitation in the present mixed spin 
chain has a gap, which ia^dfiduced via the topological 
properties of the system.ESEfP We note that for A = 1 
the phase transition to the ordered phase occurs at the 
critical value ac = 0.16. 

Next we turn to the system shown in Fig. 9. We re- 
fer to it as the "diagonal mixed-spin system" since the 



to an assembly of isolated mixed-spin chains which have 
the disordered ground state. For A = 1 the phase tran- 
sition to the ordered phase occurs at the critical value 
ac = 0.21. 

Carefully observing Fig. 10 one notices that two criti- 
cal lines intersect each other in the vicinity of a = 0.87. 
In the region a < 0.87, the area of the disordered phase 
for the diagonal system is larger than that for the colum- 
nar one. This implies that when we increase A or a in 
this region, the antiferromagnetic correlation among spin 
clusters grows up more easily in the columnar system 
than in the diagonal one, because larger spins (s = 1) 
are directly coupled with each other in the columnar case, 
and stabilize the long-range order more effectively. 

On the other hand, when a > 0.87, the situation is 
reversed, namely, the disordered phase in the columnar 
system is more stable than that for the diagonal system. 
This may be understood by observing the behavior in 
the limit a —>■ oo. In this limit— the columnar system is 
reduced to the three-leg ladder |13) while the diagonal one 
still forms the 2D network. Recall here that this three- 
leg ladder consists of two s = 1/2 chains audime s = 1 
chain, which is known to have the spin gapEa'H There- 
fore, in the limit of a ^ oo the 2D antiferromagnetic 



Quantum Phase Transitions in Two-Dimensional Spin Systems 



correlation vanishes for the columnar system, while it 
may still survive for the diagonal one. This may explain 
the behavior observed in the region a > 0.87. 

Up to now, the mixed-spin chains found experimen- 
tally _axe known to exhibit the ferrimagnetic ground 
statcliJ It may be expected that the mixed-spin chains 
with singlet ground state, as discussed here, may be also 
synthesized experimentally in the future. It would thus 
be an interesting subject to observe how such disordered 
systems are driven to the magnetically ordered phase in 
the presence of the interchain couplings, impurities, etc. 

§4. Summary 

We have investigated the quantum phase transitions in 
the 2D spin systems with ladder, plaquette, and mixed- 
spin structures. In order to quantitatively study the 
phase transitions, we have employed the systematic clus- 
ter expansion methods. It has turned out that the 
present approach improves to large extent our previous 
results on the phase diagram obtained by the non-linear 
a model as well as the modified spin wave analysis. For 
example, from the results on the ladder-structure sys- 
tems we have confirmed that the dimer expansion analy- 
sis is comparable to the QMC simulation in some param- 
eter regions. We have also studied the phase diagrams for 
the spin systems with the plaquette and the mixed-spin 
structures. For plaquette systems, starting with isolated 
plaquettes we have introduced the couplings among the 
plaquettes in two distinct ways: the one is naturally ex- 
trapolated to the 2D square lattice, while the other is 
to the 1/5 depleted square lattice. In particular, for the 
latter case, we have clarified how the plaquette chains 
studied previously are coupled to form the 1/5 depleted 
square lattice system. Also, for the mixed spin systems, 
we have considered two ways for stacking the spin chains, 
which are referred to as the columnar and diagonal sys- 
tems. It has been pointed out that the stability of the 
disordered phase non-trivially depends on how to stack 
the mixed spin chains. 

In this paper we have concentrated on the staggered 
susceptibility to establish the phase diagram. Not only 
to confirm our present results but also to get further 
information, it is desirable to calculate the elementary 
excitation spectrum in the same framework of the clus- 
ter expansion. In this direction, we have performed a 
preliminary calculation of the spin excitation gap for the 
columnar mixed spin system up to the fifth order in the 
coupling constant. The phase boundary determined in 
this way turns out to be in fairly good agreement with 
the one shown in Fig. 10. We thus believe that the 
present analysis of the staggered susceptibility, although 
it has been restricted to the fourth order, already cap- 
tures essential properties of the spin systems discussed 
in this paper. 
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